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Abstract 

Exact and asymptotic formulae are displayed for the coefficients A n used in Li's criterion for 
the Riemann Hypothesis. In particular, we argue that if (and only if) the Hypothesis is true, 
A n ~ n(Alogn + B) for n — > oo (with explicit A > and B). The approach also holds for more 
general zeta or L-functions. 

Alternative expressions are presented for the coefficients \ n introduced by Li [T] 
to recast the Riemann Hypothesis as A n > (Vn): a first representation © as a finite 
oscillatory sum, then a closely related integral representation (J2J), upon which the saddle- 
point method finally yields definite n —* +oo asymptotic estimates in one of two sharply 
distinct forms, (jHJ) or (fTTljl . depending on the falsity or truth of the Riemann Hypothesis. 
(Only the main ideas are indicated here.) 
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1 Background. 



The coefficients \ n are denned as A n = [1 — (1 — l/p) n ], or equivalently via the generating 

p 

function 



dz 



log2(— ) s^An^" 1 , (S(a) = ,(,-l)r( S /2) 7 r-/ 2 C(^)). (1) 

n=l 

The Riemann zeros are listed in pairs, as 

{p = | ± iTk}k=i,2,..., {Re T fe } positive and non- decreasing; (2) 

sums and products over zeros are performed with p and (1 — p) paired together, as usual; 
we parametrize each such pair by the single number Xk = p(l—p) = \+T\? . A "secondary" 
zeta function is 

oo 

Z{a) = ^x k - a , Re a > i. (3) 

k=l 

It extends to a meromorphic function in C, and all its poles lie at the negative half-integers 
except cr = +i p, which has the polar part [§| 

Z{\ + e) = R- 2 e- 2 + R- 1 e- 1 + 0(1)^, R-2 = ^-, R_ 1 = -^l. (4) 

2 Exact forms. 

To reexpress the X n , we use a symmetrical Hadamard product form of H(s) jU Sec. 1.10] [§], 
namely H(s) = njfeli [1 — s (l — s )/ x k], instead of the standard form chap. 12]. Hence, 

then, expanding (1 — z)~ 2 i by the generalized binomial formula, reordering in powers of 
z and identifying the output with ((TJ), we get as first result 

^-"E^Cy (6) 

Remark: other such linear relations already exist, involving the cumulants of the Stieltjes 
constants [I] and/or the sums ^ = J2 p p~i instead of Z(j); see (HI Sec. 3. 3] and references 
therein. An advantage of © is that the lesser known factors Z(j) are positive, and 
smooth (they sample the holomorphic function Z(a)). Still, this evaluation of A n involves 
cancellations that increase with n. 

An integral representation equivalent to © simply by residue calculus, but much more 
manageable, is 

K = / 1(a) ia, 1(a) = n "t,fl a - n) ZW, (7) 

7T y c r(2a + i) 

where C is a positive contour encircling (only) the subset of poles a = +1, • ■ ■ , +n of the 
integrand 1(a). 
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3 Asymptotic forms. 



The integral formula readily suggests an asymptotic (n — > oo) evaluation by the classic 
saddle-point method Sec. 2. 5]. First, the contour C is deformed in the direction of 
decreasing |/(cr)| until it crosses the nearest saddle-points of |/(ct)|. Among these saddle- 
points, the one(s) where is the largest give the dominant large-n contributions to 
the integral, through local formulae at each such saddle-point (the integrand itself may 
be asymptotically approximated as well: e.g., by a Stirling formula for T(z)). 

We caution that our subsequent assertions, involving an integrand which cannot be 
described in fully closed form, partly retain an experimental character and would warrant 
further conhrmation. 

In the present problem and for large n, the landscape of the function \I(cr)\ near the 
contour C is dominantly controlled by the T-ratio: the resulting deformation is a dilation 
of C away from the segment [l,n], then the nearest saddle-points can be of two types 
here (once n is large enough). 

1) For a on the segment (~, 1), \I(cr)\ ~ n [sin ixo T(2cr + l)]~ 1 n 2er ~ 1 Z(a) always has 
one real minimum <r r (n) (tending to ~ as n — > oo; other real saddle-points lie below a = ~, 
and are sub dominant). 

2) There may exist complex saddle-points er c (n), with imaginary parts proportional 
to n; we may focus just on the upper half-plane, as the lower half-plane gives complex- 
conjugate ("c.c") contributions. Wherever Z(a) is dominated by the term Xk~ a from one 
Riemann zero Xk (naively, xi), the saddle-point equation (in the Stirling approximation 
for the T-ratio) is = 4- log \ ~ log(cr 2 — n 2 ) — 2 log 2a — logx^, yielding the formal 
saddle-point location 



Here, any zero Xk on the critical axis (rfc real) yields a purely imaginary a c (n), not 
eligible: it lies outside the domain of convergence of (j3J), and its contribution would be 
subdominant in any case. The discussion then fundamentally splits depending on the 
presence or absence of zeros off the critical axis. 

[RH false] If the zero Xk lies off the critical axis (and selecting argTfc > 0), the as- 
sociated complex saddle-point cr c (n) is relevant: it lies inside the domain of convergence 
{Re a > |} as soon as n | Im l/r^l > 1, and it will exponentially dominate the real saddle- 
point (as seen later): by the standard calculation, its asymptotic contribution to A n is 
[{r k + i/2)/(r fc — i/2)] n , which grows exponentially in modulus and fluctuates in phase. 
This result can also be rigorously confirmed directly from (£Q): by a conformal mapping ar- 
gument [T], the function 4- logH(Y^) has precisely the points Zk = {j k — i/2)(r fc + i/2) -1 
and z_k — ~Zk as simple poles of residue 1 in the unit disk; by a general Darboux theorem 
|H1 chap. VII §2], this implies the asymptotic form J2k( z ±k)~ n f° r the Taylor coefficients 
A„ of that function, i.e., 



cr c (n) = ni/2r k . 





[ + c.c] 




n 



oo. 



[RH false] 



(9) 
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(An infinite set of these zeros poses no extra problem: their contributions form an asymp- 
totic sequence.) 

On the other hand, the Darboux approach does not resolve the case [RH true], in 
which the poles z±k all lie at the same (unit) distance, and cluster at z — 1 ! 

[RH true] By contrast, the saddle-point analysis of (J7|) still appears to work. Now all 
the Tfc are real, Z(a) = 0(Z(Re a) | Im cr|~ 3 / 2 ) in {Re a > |}, and the contour C can be 
freely moved towards the boundary of the half-plane without crossing any of the a c (n) (all 
purely imaginary). Hence the only dominant saddle-point in this case is a r (n) e 
it is shaped by the double pole of Z(a) at \ (itself generated by the totality of Riemann 
zeros), in the form a r (n) ~ | + One may then proceed as usual in the quadratic 

approximation of log 1(a) around a T (n), but this is not so fit for a confluent case (a = | 
is a singular point). Here, it is at once simpler and more accurate to keep on deforming a 
portion of the contour C nearest to a — \ until it fully encircles this pole (now clockwise), 
and to note that the resulting additions to the integral are asymptotically negligible. 
Hence for [RH true], the result is 

T(a + n)T(a -n) 



X n ~ (-l) n 2nRe S(T=1/2 



-Z(a) ~27rn[2i2_ 2 log7i-2i2_ 2 (l- 7 ) + iL_ 1 ]; 
1 (2a + 1) J 

(10) 

(7 = Euler's constant); specifically for the Riemann zeros, using the explicit values (J3J), 

A n ~ |n(logn — log27r — 1 + 7), n — > 00. [RH true] (11) 

The form (fTU|) covers more general situations: the zeros of, e.g., Dedekind zeta functions 
and some Dirichlet L-functions ^UI, i n their [RH true] case (i?_ 2 is always positive, in 
agreement with Li's criterion). 

We see a good agreement of (fTTjl with numerical data |H] for n < 3300 - still a bit 
short to give full confidence in (|11|). however. (This agreement is even improved in the 
mean if we include the contribution like (fTU|) but from the next pole, 



S\ = (-l) n 2nRes n 



T(a + n)T(a -n) 

-Z(a) 



2Z(0) = +1, (12) 



T(2a + 1) 

although this term should not be asymptotically meaningful - larger oscillatory terms are 
also present.) 

Finally, the result becomes even simpler for some pure linear combinations of the 
Stieltjes cumulants. The latter may be defined in degree n as 

9n = (- 1 ) T^ lo S S C + s Uo = - Jim <y > 13 

ds n m^+oo I ^— ' m n 

v m=l J 



(cf. formula (4.1) in [1: this relates to the Euler factorization of ((s) over the primes). 
Our (y£ corresponds to (— l) n (n — 1) ! 77^—1 in |Q]: e.g., g\ = 7 = — T] . Theorem 2 in P 

evaluates the differences (A n — S n ) where S n = — ^2 (™) r]j_i = ^2 (j-iy. (j) 9j- ^ s a 

consequence for n — >• 00, (fTTj) implies 

5 n = o(n); [RH true] (14) 
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whereas in the opposite case, we expect S n ~ A n as in (jHJ): it will oscillate between 
exponentially large values, negative and positive, but appreciable in absolute size only 
beyond n pa min {[Im l/r^] -1 }. 

{argr fc >0} 
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